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Abstract 

A measurement is presented of inelastic photo- and electroproduction of J/ tp mesons in ep 
scattering at HERA. The data were recorded with the HI detector in the period from 2004 to 
2007. Single and double differential cross sections are determined and the helicity distribu- 
tions of the J/ip mesons are analysed. The results are compared to theoretical predictions 
in the colour singlet model and in the framework of non-relativistic QCD. Calculations in 
the colour singlet model using a kx factorisation ansatz are able to give a good description 
of the data, while colour singlet model calculations to next-to-leading order in collinear 
factorisation underestimate the data. 
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1 Introduction 



The description of the process of charmonium production in interactions of photons and hadrons 
is a challenge to theory, since it involves both the production of the heavy quark system and the 
formation of the bound state. Charmonium production in electrorQ-proton collisions at HERA is 
dominated by photon-gluon fusion: a photon emitted from the incoming electron interacts with 
a gluon from the proton to produce a cc pair that evolves into a charmonium state. In the colour 
singlet model, only those states with the same quantum numbers as the resulting charmonium 
contribute to the formation of a bound cc state. This is achieved by radiating a hard gluon in a 
perturbative process. In the factorisation ansatz of non-relativistic quantum chromodynamics, 
also colour octet cc states contribute to the charmonium production cross section via soft gluon 
radiation. 

Previous measurements in electroproduction (ep) and photoproduction (^yp) at HERA 03-0 are 
not described by predictions in the colour singlet model to leading order. In contrast, the calcu- 
lation of photoproduction cross sections to next-to-leading order (NLO) (H showed a reasonable 
description of the photoproduction cross sections. The calculation proved that the corrections 
with respect to leading order results are very large, increasing towards large transverse mo- 
mentum of the J/tp meson. The same calculation, repeated recently with an up-to-date set of 
theoretical parameters flU, results in a prediction which is about a factor of three below the 
measured cross sections, indicating that corrections beyond NLO are needed and/or that contri- 
butions from colour octet states may be sizable. 

In this paper a measurement is presented of inelastic J/tp meson production at HERA. The 
measurement uses a larger data sample than previous results 03-|4l and benefits from improved 
systematics. The data sets were collected in the years 2004 to 2007 with the HI detector. The 
J/tp meson candidates are identified by the leptonic decay into two muons or electrons. The 
cross sections are measured for both electroproduction and photoproduction. For the photo- 
production sample J/ip meson polarisation variables are determined. The data samples are 
restricted to the region of phase space where contributions from diffractive charmonium pro- 
duction are suppressed. 



2 Theoretical Models 

In order to describe inelastic charmonium production in the framework of perturbative QCD 
different models have been proposed, such as the colour-evaporation model H101I11] . the colour- 
singlet model (CSM) [fT2]4T6l . the factorisation ansatz in non-relativistic quantum chromody- 
namics (NRQCD) [fT7]4T9l and soft colour interactions [|20l . In this paper the most recent cal- 
culations using the CSM or NRQCD are compared to the data. 

In the CSM, only charm quark pairs in a colour singlet state with the same quantum numbers 
as the resulting charmonium contribute to the formation of a bound cc state. This is achieved 
by radiating a hard gluon in the perturbative process. The factorisation ansatz in NRQCD 
includes also colour octet cc states in the charmonium production cross section. The size of 

'In this paper "electron" is used to denote both electron and positron. 
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these colour octet contributions, described by long distance matrix elements (LDME), is defined 
by additional free parameters which were determined in fits to the Tevatron data [12111 . The 
NRQCD factorisation approach contains also the colour singlet model which is recovered in the 
limit in which the colour-octet LDME tend to zero. 

The following calculations are compared to the measurements presented in this paper: 

• A calculation of J ftp meson photoproduction via a colour singlet mechanism [9] provides 
predictions for both cross sections and helicity distributions to next-to-leading order. The 
uncertainty of this calculation is estimated by variations of the charm quark mass and the 
factorisation and renormalisation scales. 

• A calculation at NLO for photoproduction cross sections includes the full framework of 
NRQCD E2l . The uncertainty of this calculation is dominated by the limited knowledge 
of the LDMEs. 

• CSM predictions in the kr factorisation approach are employed as implemented in the 
MC generator CASCADE [|23l . Higher order parton emissions based on the CCFM evolu- 
tion equations ||24l| are matched to 0(a s ) matrix elements in which the incoming parton 
can be off-shell. The uncertainty on the calculation is estimated by varying the renormal- 
isation scale by a factor of two. In addition polarisation variables in the k T factorisation 
approach are calculated analytically E51 . 

Parameters and variations used in the theoretical calculations are given in tabled] 

3 HI Detector 

The HI detector is described in detail elsewhere ll28l . Here only the components essential to the 
present analysis are briefly described. A right handed Cartesian coordinate system is used with 
the origin at the nominal primary ep interaction vertex. The proton beam direction defines the 
z axis. The polar angles 9 and transverse momenta P T of all particles are defined with respect 
to this axis. The azimuthal angle defines the particle direction in the transverse plane. The 
pseudorapidity is defined as r\ = — In tan |. 

Charged particles emerging from the ep interaction region are measured by the central track- 
ing detector (CTD) in the pseudo-rapidity range \q\ < 1.74. The CTD consists of two large 
cylindrical central jet drift chambers (CJC) which are interleaved by a 2-chamber and arranged 
concentrically around the beam-line in a magnetic field of 1.16 T. The CTD provides triggering 
information based on track segments from the CJC Il29ll30l . and on the ^-position of the vertex 
from the 5-layer multi-wire proportional chamber Oil which is situated inside the inner CJC. 
To provide the best possible spatial track reconstruction, CTD tracks are linked to hits in the 
vertex detector, the central silicon tracker CST 11321 . The CST is installed close to the interac- 
tion point, surrounding the beam pipe in the pseudo-rapidity range \q\ < 1.3 and consists of 
two layers of double sided silicon strip sensors. 

Charged and neutral particles are measured in the liquid argon calorimeter (LAr) ll33l which sur- 
rounds the tracking chambers and covers the range — 1.5 < r\ < 3.4 and a lead/scintillating-fibre 
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CSM (k T factorisation), CASCADE 11231 
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PDF 
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A (3) 
iv QCD 



CCFM set AO [|27l 

Ho = 



4 



mfi. + P 



2 

T,ip 



1.5 GeV 
200 MeV 



Table 1: Summary of the parameters employed in the CSM and NRQCD calculations used 
to compare to the measurements in this paper. In this table PDF means parton distribution 
function of the proton, s denotes the invariant mass square of the hard subprocess and Q± the 
initial transverse momentum of the partonic system (7g). 

calorimeter SpaCal [|34l . covering the backward region —4.0 < r] < —1.4. The calorimeters 
are surrounded by the solenoidal magnet and the iron return yoke. The yoke is instrumented 
with 16 layers of limited streamer tubes, forming the central muon detector (CMD) in the range 
-2.5 < r] < 3.4. 

The luminosity determination is based on the measurement of the Bethe-Heitler process ep — > 
ep7, where the photon is detected in a calorimeter located downstream of the interaction point 
in the electron beam direction at z = —104 m. 
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Figure 1 : Invariant mass spectra of two oppositely charged leptons after all selection cuts for a) 
the photoproduction (jp) sample and b,c) the electroproduction (ep) samples as listed in table[2] 
The lines show the results of fits to signals and backgrounds. 



4 Data Analysis 

The kinematics of inelastic charmonium production at HERA are described using the following 
variables: the square of the ep centre of mass energy s = (p + k) 2 , where p and k denote the 
four vectors of electron and proton respectively; the negative squared four momentum transfer 
Q 2 = —q 2 , where q is the four vector of the virtual photon; and the mass of the hadronic final 
state W rp = a/G» + q) 2 . W 1P is related to the scaled energy transfer y = (p ■ q)/(p ■ k) via 
Wryp 2 = ys — Q 2 . In addition, the elasticity of the J/ij) meson production process is defined as 
z = (pip • p)/(q • p), where p^ is the four momentum of the J/ip meson. The elasticity denotes 
the fractional energy of the photon transferred to the J/ij) meson in the proton rest system. 

Events are selected separately in the photoproduction and electroproduction regimes. Photo- 
production events are selected by requiring that no isolated high energy electromagnetic cluster, 
consistent with a signal from a scattered electron, is detected in the calorimeters. This limits the 
virtuality to values of Q 2 < 2.5 GeV 2 , resulting in a mean value of (Q 2 ) 0.085 GeV 2 . Con- 
versely, for the electroproduction sample, a scattered electron with energy of more than 10 GeV 
is required to be reconstructed in the backward calorimeter (SpaCal), corresponding to a range 
in photon virtuality 3.6 < Q 2 < 100 GeV 2 . 

In this analysis the photon virtuality Q 2 is reconstructed from the scattered electron energy 
E' e and polar angle 9' e as Q 2 = AE e E' e cos 2 (9' e /2), where E e denotes the energy of the beam 
electron. The variable y is reconstructed using the relation y = ^2 h (E — p z )/2E e for pho- 
toproduction lf35l and y = J2h(E — Pz)/z~2(^ — Pz) f° r electroproduction [|36ll . The sums 
in the numerator include all particles of the hadronic final state without the scattered electron, 
which is only included in the sum of the denominator for electroproduction. The elasticity z is 
then obtained from z — (E — p z )j/^p/ ^ h {E — p z ), where (E — p z )j/^p is calculated from the 
decay particles of the J/ij) meson. The kinematics of the final state particles are obtained from 
charged particle tracks reconstructed in the CTD and energy depositions in the LAr and SpaCal 
calorimeters H37H38H . 

The J/ip meson candidates are reconstructed through their decays into two oppositely charged 
muons or electrons. These decay leptons are reconstructed as charged particles in the CTD 
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Photoproduction Electroproduction 




J /if) — > J /if) — > n + ji~ J /if) — > e + e~ 


kinematic range 




Q 2 < 2.5 GeV 2 3.6 < Q 2 < 100 GeV 2 




60 GeV < W JP < 240 GeV 




P T ^ > 1 GeV P* 4 , > 1 GeV 




0.3 < z < 0.9 


event selection 




P T/ > 800 MeV 




20° < Bp < 160° 20° <9 fl < 160° 20° < 9 e < 150° 




Ntyi, > 5 (in the range 20° < 9 < 160°) 


event samples 




2320 ±54 501 ±34 290 ± 24 


r 

'"int 


165 pb- 1 315 pb- 1 315 pb- 1 



Table 2: List of selection cuts and event yields for each of the three data samples. 

with a transverse momentum of at least 800 MeV. Muon candidates are identified as mini- 
mum ionising particles in the LAr calorimeter or through track segments in the CMD (20° < 
9^ < 160°) [|39l . Electron candidates are identified through their energy deposit in the central 
calorimeter (20° < 9 e < 150°) Q40] - For trigger reasons the photoproduction sample of J/if> 
meson events is restricted to decays into while the electroproduction sample includes 

both leptonic decay channels. The photoproduction sample was recorded in the years 2006 and 
2007 and corresponds to an integrated luminosity of C = 165 pb -1 , while the electroproduction 
sample was recorded in the years 2004 to 2007 and corresponds to an integrated luminosity of 
£ = SISpb- 1 . 

The measurement is performed in the kinematic range 60 < W 1V < 240 GeV, 0.3 < z < 0.9 
and Pt,i/}(Pt i/>) > 1 GeV. In photoproduction the transverse momentum Pt,^> is measured in 
the lab frame, while in electroproduction the transverse momentum P£ . is calculated in the Yp 
rest frame. To suppress contributions from diffractive production of J /if) and ip(2S) mesons, 
selected events are required to contain at least five reconstructed tracks in the central region of 
the detector (20° < 9 < 160°). The reconstruction efficiency accounts for this experimental cut 
and the measured cross sections are corrected for this track multiplicity cut. 

Figure \T\ shows the invariant mass spectra of the leptons in the selected event samples. The 
number of signal events, Nj/^, is obtained in all bins of the cross section measurements from 
a fit to the mass distributions in the interval 2 < mu < 6 GeV. For the decay into muons 
the signal peak is described using a modified Gaussian [7|. In the case of a decay into two 
electrons an exponential is added to the lower mass region of the signal Gaussian in order to 
account for the radiative tail [1391 . For the differential cross section measurements, the width 
and asymmetry term of the mass peak in each bin are fixed to the values obtained from the full 
samples. For both decay channels, the background is parametrised by a polynomial of third 
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order. At m u 3.7GeV, the nominal mass of ip(2S) mesons, an additional Gaussian with 
fixed position and width is allowed in all analysis bins. 

The selection criteria and the obtained event samples are summarised in tabled 



5 Monte Carlo Simulations 

Cross sections and polarisation parameters are derived by correcting the measured number of 
events and angular distributions for detector effects, such as detector resolutions and ineffi- 
ciencies. Several Monte Carlo generator programs are used to determine the corrections. All 
samples are passed through a detailed simulation of the HI detector response based on the 
GEANT program ll43l and through the same reconstruction and analysis algorithms as used for 
the data. 

Signal events are generated using the Monte Carlo generator CASCADE [|23l . Elastic and 
proton-dissociative production of ip(2S) mesons is simulated using DIFFVM [431 with pa- 
rameters tuned to describe the results of previous HI measurements Il45ll46ll . The Monte Carlo 
generator Pythia ll47l is used for the description of the contribution from b hadron decays as 
described in section |6] All generators use the JETSET part of the Pythia program ll47ll to 
simulate the hadronisation and decay processes. 

Signal events as simulated with the Monte Carlo generator CASCADE are compared with the 
data after final selection in figures |2] and [3j All data distributions in these figures are cor- 
rected for contributions from non-resonant background events using a sideband method de- 
scribed in ||39l . 

Corrections as a function of W 1V and Pt^ in bins of the elasticity z are applied to the CASCADE 
Monte Carlo simulation in order to describe the data. Details of the procedure are described 
in [|39l . In figure |2] distributions for the photoproduction sample are compared to CASCADE 
Monte Carlo predictions before and after correction for the observables Pt,^ Pr,ip, #v ^yp 
and z. Similarly, in figure |3l the summed distributions for the two electroproduction samples 
(J ftp — > and J/ip — > e + e~) are shown for the observables P£ 9^, Q 2 , X-P^charged, W 1P 

and z. Here, £-Pr,charged is the scalar sum over the transverse momenta of all measured charged 
particles except for the scattered electron and the J/ip meson decay leptons. The corrected 
Cascade simulation gives a good description of all aspects of the data and is used to correct 
the data for losses due to limited acceptance and efficiency of the detector. 



6 Backgrounds 

Remaining backgrounds to prompt J/ip meson production in the selected sample originate from 
feed-down processes, i.e. J/ip mesons produced in decays of diffractively or inelastically pro- 
duced ip(2S) mesons and \c mesons or of b hadrons. 

Inelastic production of ip(2S) mesons with a subsequent decay into J/ip mesons is expected to 
contribute about 15 — 20 % to the selected J/ip meson samples (USSl. Since the production 
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processes are the same, the inelastic ip{2S) mesons show similar dependences on the kinematic 
variables. 

Diffractive production of if)(2S) mesons contributes at large values of z by decays into a J/tp 
meson and two charged pions. These events typically contain three or four reconstructed 
charged tracks in the central detector (20° < 9 < 160°). In figure Hk) the distribution of 
the charged track multiplicity measured in the central detector is shown for the photoproduction 
sample selected using all selection criteria given in table [2] except for the track multiplicity cut, 
which is relaxed to N Trk > 3. The data are described by the sum of the CASCADE simulation 
and the prediction for diffractive ip(2S) production, as simulated using the DIFFVM Monte 
Carlo generator. In the final selection remaining contributions from diffractive ip(2S) meson 
production amount to about 1.3% in the total sample and to about 5% in the highest elasticity 
bin, 0.75 < z < 0.9. 

The fraction of events arising from b hadrons decaying into J/ip + X is estimated using the 
PYTHIA simulation. The Pythia prediction is scaled by a factor of 2, based on results from pre- 
vious measurements of beauty production at HERA H49U50I . This scaled prediction by Pythia 
amounts to 5% in the total sample and about 20% in the lowest z bin. It is confirmed within 
uncertainties by the following determination using data. The fraction of events in the photo- 
production sample containing b hadrons is estimated using the impact parameter of the decay 
muons to exploit the lifetime signature of b hadrons. The impact parameter, 5, of the decay 
muon tracks is defined as the distance of closest approach in the transverse plane to the recon- 
structed primary vertex. The sign of the impact parameter is defined as positive if the angle 
between the decay muon and the J ftp meson momentum direction is less than 90 °, and is de- 
fined as negative otherwise. A signed significance S = 5/a(S) is reconstructed by weighting 
the reconstructed signed impact parameter with its uncertainty I15T1 . Figure lib) shows the dis- 
tribution of the signed significance for events in the interval 0.3 < z < 0.4. The histogram is 
filled with the signed significance of the decay muons for all events where both muon tracks 
have at least one hit in the CST The fraction of events coming from the decay of b hadrons is 
obtained from a fit of the significance distribution of CASCADE (simulating prompt J/ip meson 
production) plus Pythia (simulating bb events with subsequent decays into J/ip + X) to that of 
the data. The fit results are dominated by the region of small signed significances, S < 3, due to 
large statistical uncertainties at larger values of S. The distribution of the data is corrected for 
non-resonant contributions using the side bands |f3~9ll . The relative contribution from b hadrons 
as resulting from the fits are shown in figure St) for three bins of z. The scaled predictions 
from Pythia are found to be in good agreement with the measured fractions, indicating that 
the background from b hadrons is under control. 

The contribution from \c production and decay was studied [3] and found negligibly small in 
the present kinematic region, 0.3 < z < 0.9. 



7 Systematic Uncertainties 

The sources of systematic uncertainties of the cross section measurement are listed in table [3] 
and are detailed in the following: 
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Source 




Uncertainty [%] 






Photoproduction Electroproduction 




J /if) —> 


J /if) —> fi + fi 


J /if) -> e + e~ 


Decay leptons reconstruction 


1 


1 


2 


Decay leptons identification 


Q 
o 


o 


Q 
o 


Number of signal events 


2 


2 


4 


Trigger 


3 


2 


2 


Scattered electron energy scale 




2 


2 


Hadronic final state energy scale 


4 


3 


3 


Integrated luminosity 


4 


3.2 


3.2 


Model uncertainties 


5 


5 


5 


Decay branching ratio 


1 


1 


1 


Sum 


9.0 


8.2 


9.1 



Table 3: Systematic uncertainties of the J /if) meson production cross section. The total system- 
atic uncertainty is the sum of the contributions added in quadrature. 



• The uncertainty on the cross section due to the track and vertex reconstruction efficiency 
has been determined to be 1 % for J /if) — > /i/i and 2 % for J /if) — > ee. 

• The efficiency for the identification of the leptons is determined using a high statistics 
sample of events of elastically produced J /if) mesons [3~9l . The detector simulation is 
reweighted to match the efficiency measured in the data as necessary. Remaining differ- 
ences are smaller than 3 % everywhere and are taken as systematic uncertainty. 

• The systematic uncertainty on the determination of the number of signal events, obtained 
by a fit to the mass distributions in every analysis bin, is determined by a variation of the 
extraction method. Comparing the number of signal events for binned and unbinned log- 
likelihood fits yields a systematic uncertainty of 0.5 %. In addition, the result from the 
fit to background and signal is compared to the number of signal events above the fitted 
background function in the mass window between 2.95 and 3.2 GeV. An uncertainty of 
2 % for the decay into muons and 4 % for the electrons is found. The uncertainty for the 
electron is larger due to an additional uncertainty originating from the description of the 
radiative tail. 

• The trigger efficiencies are determined using independent trigger channels. For the elec- 
troproduction sample the trigger efficiency is measured to be (97 ± 2) %. In the pho- 
toproduction sample the trigger efficiency depends mainly on the identification of the 
decay muons in the central muon system. The efficiency amounts to about 70 % with a 
systematic uncertainty of 3 %. A detailed description of the determination of the trigger 
efficiencies can be found in 11391 . 

• For the electroproduction sample the measurement of the scattered electron energy is 
known with a scale uncertainty of 1 %. The uncertainty of the scattering angle is 1 mrad. 
Both uncertainties combined lead to an uncertainty of the cross section measurement of 
2 % on average. 
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• The hadronic energy scale uncertainty is 4 % in the LAr and 7 % in the SpaCal. This 
leads to an uncertainty on the cross sections measurement of 3 % for the electroproduction 
sample and 4 % for the photoproduction sample. 

• The integrated luminosity is known to a precision of 3.2 % for the electroproduction sam- 
ple and 4.0 % for the photoproduction sample. 

• The dependence of the result on model assumptions made in the CASCADE Monte Carlo 
simulation were investigated and found to amount to 5 % in total. The model uncertainty 
arising from the knowledge of the decay angular distributions, explained in section |9l is 
determined by variation of the parameter a in the simulation by ±0.3. This variation 
results in a change of the cross section of up to 4 %. The systematic uncertainty originat- 
ing from the uncertainty of the slope of the P Ty ^ (P£ ^) distribution in the simulation is 
determined by a variation of this distribution as described in [|39ll . This variation results 
in a change of the cross section of up to 4 %. 

• The branching ratios of the leptonic decay channels of the J ftp meson are known with an 
accuracy of 1% 03- 

The total systematic uncertainty is obtained by adding all the above contributions in quadrature. 
A total systematic uncertainty of 9 % is determined for the photoproduction sample. For the 
combined electroproduction cross section the total systematic uncertainty is 8.5% The same 
uncertainties are attributed to all bins of the cross section measurement. For the measurement of 
the helicity distributions only the uncorrected systematic uncertainties are taken into account. 
They amount to about 3.5 % and are negligible compared to the statistical uncertainties. 



8 Cross Section Measurements 

The cross section measurement is performed in the kinematic range 60 < W 7P < 240 GeV, 
0.3 < z < 0.9 and Pt,ip(Pt,^) > 1 GeV. The photon virtuality Q 2 is limited in the electro- 
production analysis to 3.6 < Q 2 < 100 GeV 2 and for the photoproduction sample to Q 2 < 

2.5 GeV 2 . 

For the measurement of differential cross sections the number of signal events in each bin is 
corrected for detector inefficiencies and acceptance and normalised to integrated luminosity and 
branching ratio. They are not corrected for QED radiative effects. The electroproduction cross 
sections, measured from J ftp — > fifi and J/ip — » ee, are combined 11391 . The differential cross 
sections are bin-centre corrected using MC simulations. In order to avoid model dependencies, 
the measured cross sections are not corrected for contributions from backgrounds as described 
in section[6] All measured cross sections are listed in tables 0]-[9] together with statistical and 
systematic uncertainties. 

For the photoproduction sample the measured ep cross sections are transformed to 7p cross sec- 
tions using the photon flux factors presented in table [51 calculated in the Weizsacker Williams 
approximation [|52l . The differential J/ip meson photoproduction cross section is measured 
as function of the elasticity z and the squared transverse momentum P| , of the J/ip meson. 
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The total ^p cross section is measured in bins of the photon proton centre of mass energy 
Wryp. The results are displayed in figure [5]and show a reasonable agreement with the prediction 
from the CASCADE MC generator. A variation of the renormalisation scale by a factor of two 
(0.5/xo < A*r < 2/i ) has little effect as shown by the band in the figures. In addition to the 
CASCADE prediction, the remaining contributions from diffractive ip(2S) mesons and from b 
hadrons are shown. The distributions in P| ^ and z are further investigated by dividing the sam- 
ple into bins of P^ and z, respectively as shown in figure |6j The z distribution tends to flatten 
off towards larger values of Pt ^ presented in figure [6^). It can be seen that differences between 
the data and the CASCADE prediction are localised at low elasticities and low transverse mo- 
menta of the J/ip mesons, where CASCADE overshoots the data, and at large elasticities and 
large transverse momenta, where CASCADE is below the data. Taking into account that the 
measured cross section in the lowest elasticity bin includes a significant fraction of about 20 % 
of events originating from b hadron decays, the difference to the CASCADE prediction is even 
more significant. 

Results for electroproduction are shown in figure [7] and figure [8] Differential ep cross sections 
are measured as functions of the photon virtuality Q 2 , the squared transverse momentum of 
the J ftp meson in the photon proton rest frame P^, the energy W 1V and the elasticity z. 
Figure [8] shows differential cross sections as a function of the elasticity z in bins of and 
as a function of P^ in bins of z. A comparison of the electroproduction data with predictions 
from the Monte Carlo generator CASCADE reveals in general a reasonable agreement with the 
data. Differences in shape can be seen in the differential cross section as a function of P^%- 

For photoproduction, several theory calculations to next-to-leading order have been performed 
and are compared with the data in figure [9] A calculation in the CSM at NLO [8 J was repeated 
using up-to-date sets of scale parameters J9l|22l, yielding predictions as shown in figure Hk)-b). 
The shapes of the data are reasonably described, whereas the normalisation of the prediction is 
about a factor three below the data, with large uncertainties, indicating that corrections beyond 
next-to-leading order are necessary in order to describe the data. Estimates of the NNLO con- 
tribution for charmonium production at the Tevatron 115311541 indicate that these contributions 
can be large indeed. 

The calculation to next-to-leading order has been extended to include colour octet contributions 
resulting in a larger cross section [22J. A comparison of this prediction with the data is shown in 
figure[9]:)-d). The dominant uncertainty arises from the difference in the predicted cross section 
when using LO colour octet LDMEs or higher order improved LDMEs [|22l . The NRQCD 
prediction fails however in describing the shape of the differential cross section as a function of 
the elasticity z, even within the presently large uncertainties of the calculation. 

9 Polarisation Measurement 

The measurement of the J/ip meson helicity distributions provides an independent method to 
distinguish between different production mechanisms. The measurement is performed for the 
photoproduction data sample. The J ftp meson polarisation is measured by analysing the decay 
angle distributions of the J/ip meson, and their dependence on Pt^ and z, in two comple- 
mentary frames ll55l : the helicity frame and the Collins-Soper frame. In the helicity frame the 
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polarisation axis z in the J/ip meson rest frame is defined by the flight direction of the J/ip 
meson in the 7^ rest frame, whereas the polarisation in the Collins-Soper frame is measured 
with respect to the bisector of proton (— p p ) and photon (jiy) in the J/ip meson rest frame Il56ll . 
Subsequently, the frame-dependent polarisation axis is taken as z axis of a right handed coordi- 
nate system, where the x and z axis lie in a plane spanned by the photon and proton directions. 
The y axis is perpendicular to this plane and is the same in both reference frames. The polar 
(9*) and azimuthal (</>*) angles of the positive decay muons are used. 

The parametrisation of the measured decay angle distributions as function of cos(#*) and (p* is 
given by (EH: 



da 



oc l + acos 2 6>*; (1) 



d cos 9* 



da a v 

w cx l + j + 5 »V. CD 

The polarisation variables a and v can be related to elements of the spin density matrix for 
the J/ip meson. Moreover, a = +1 and —1 corresponds to fully transverse and longitudinal 
polarisation of the J/ip meson, respectively. 

A x 2 At is performed in each bin of the polarisation measurement, comparing data to Monte 
Carlo samples on reconstruction level probing values for a and v between —1 and +1. Sys- 
tematic uncertainties on this measurement are negligible compared to rather large statistical 
uncertainties. The results for a and v as a function of P T ^ and z are presented for the helicity 
frame in figure [TT] and in figure [[2] for the Collins-Soper frame. The values for the polarisation 
parameters in both frames are listed in table [TOl 

Within uncertainties the J/ip mesons produced inelastically at HERA are unpolarised. The mea- 
surements are compared to predictions using a kx factorisation ansatz [25J and to calculations in 
the CSM in collinear factorisation at leading order E51 and next-to-leading order [9]. The pre- 
dictions in the k T factorisation ansatz describe the data. The NLO calculations show a similar 
trend within large uncertainties. In contrast, the leading order CSM calculation predicts larger 
values for the polarisation variables than the measured ones for many bins and is disfavoured 
by the measurement. A similar measurement was published by the ZEUS collaboration in a 
different kinematic range 11571 . 



10 Conclusions 

A measurement of inelastic J/ip meson production is performed. Differential cross sections 
with improved statistical and systematic uncertainties are presented for both electroproduction 
and photoproduction. Polarisation parameters for the photoproduction of J/ip mesons are mea- 
sured in two different reference frames, the helicity frame and the Collins-Soper frame. 

The data are compared to a number of recent theory predictions. It is found that predictions 
based on k T factorisation in the colour singlet model are able to describe the cross sections and 
the helicity distributions well. Calculations based on collinear factorisation in the colour singlet 
model at next-to-leading order produce a reasonable description of the shape of the measured 
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cross sections, but are lower in normalisation. They give an acceptable description of the polar- 
isation parameter measurements within the large uncertainties. The failure to describe the cross 
section measurements and the strong sensitivity to scale variations indicate that calculations be- 
yond next-to-leading order are necessary. Moreover contributions from colour octet states may 
be significant. 
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Figure 2: Control distributions of the photoproduction sample: a) the transverse momentum 
Pt^ of the muon tracks, b) the polar angle # M of the muon tracks, c) the transverse momentum 
Pt,iP of the J ftp meson, d) the polar angle 9^ of the J ftp meson, e) the elasticity z and f) the 
photon proton centre-of-mass energy W lp . The data are compared with predictions from the 
corrected CASCADE Monte Carlo simulation (solid lines), normalised to the number of entries 
in the data. The uncorrected CASCADE Monte Carlo prediction is shown as dashed line. 
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Figure 3: Control distributions of the electroproduction sample: a) The squared transverse 
momentum of the J/ij) meson in the photon proton rest frame P£%, b) the polar angle of the 
J/ip meson 6^, c) the photon virtuality Q 2 , d) the scalar transverse sum EP^charged, e) the 
elasticity z and f) the photon proton centre of mass energy W 7P . The data are compared with 
predictions from the corrected CASCADE Monte Carlo simulation (solid lines), normalised to 
the number of entries in the data. The uncorrected CASCADE Monte Carlo prediction is shown 
as dashed line. 



21 



Inelastic J/ip Photoproduction 




0.1 0.2 0.3 0.4 



Figure 4: a) Distribution of the multiplicity of tracks, Nxrk, m the central region of the detector 
(20° < 9 < 160°) for the photoproduction sample, b) signed significance distribution S for the 
photoproduction sample at low elasticities (0.3 < z < 0.4) and c) measured contribution from 
b hadron decays for three bins of the elasticity z in comparison with the prediction based on 
Pythia (scaled up by a factor of two) and CASCADE. 
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Figure 5: Differential meson photoproduction cross sections for the kinematic range 60 < 
W 7P < 240 GeV, 0.3 < z < 0.9 and P T ^ > 1 GeV, as functions of a) the elasticity z, 
b) the photon proton centre of mass energy W 1P and c) the squared transverse momentum of 
the J/ip meson P?^- The inner error bar represents the statistical uncertainty and the outer 
error bar indicates the statistical and systematic uncertainties added in quadrature. The data are 
compared to the predictions from CASCADE (solid line). The uncertainty band of the CASCADE 
prediction arises from a scale variation by a factor of two. The dashed and dotted lines indicate 
the remaining background from diffractive ift(2S) or b hadron decays respectively as estimated 
using MC simulations. 
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Inelastic J/i/j Photoproduction 
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Figure 6: a) Differential J/ip meson cross sections as a function of z in four bins of Pt^ and 
b) differential J/ip meson cross sections as a function of in four bins of z. The inner 
error bar represents the statistical uncertainty and the outer error bar indicates the statistical and 
systematic uncertainties added in quadrature. For visibility, the measured cross sections are 
scaled by the factors indicated in the figures. The data are compared to the predictions from 
Cascade (lines). 
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Figure 7: Differential J / ip meson cross sections for the kinematic range 3.6 < Q 2 < 100 GeV 2 , 
60 < W 1P < 240 GeV, 0.3 < z < 0.9 and P*^ > 1 GeV, as functions of a) the photon 
virtuality Q 2 , b) the squared transverse momentum of the J ftp meson in the photon proton rest 
frame P^, c) the energy in the photon proton rest frame W 1V and d) the elasticity z. The inner 
error bar represents the statistical uncertainty and the outer error bar indicates the statistical and 
systematic uncertainties added in quadrature. The data are compared to the predictions from 
Cascade (solid line). The uncertainty band of the Cascade prediction arises from a scale 
variation by a factor of two. 
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Figure 8: a) Differential J/^ meson cross sections as a function of z in three bins of P^ and 
b) differential cross section as a function of P^ in three bins of z. For visibility, the measured 
cross sections are scaled by the factors indicated in the figures. Predictions from CASCADE are 
shown as solid line. 
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Inelastic J/ip Photoproduction 

a) b) 




Figure 9: Differential J/ip meson photoproduction cross sections for the kinematic range 60 < 
W JP < 240 GeV, 0.3 < z < 0.9 and Pt^ > 1 GeV as functions of the squared transverse 
momentum of the J/ip meson Pj,, (a) and c)) and the elasticity z (b) and d)). The inner 
error bar represents the statistical uncertainty and the outer error bar indicates the statistical 
and systematic uncertainties added in quadrature. The data are compared with calculations to 
next-to-leading order: a,b) a colour singlet model (CSM) calculation and c,d) a NRQCD 
calculation including contributions from colour octet states (CS + CO) ll2~2~ll . The colour singlet 
component (CS) of the latter calculation is shown separately in addition. 
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Figure 10: Distributions for a) cos(0*) in the range 0.6 < z < 0.75 and b) (f>* in the range 
1 < Ptj, < 2 GeV. The data are compared to the corrected CASCADE simulation with three 
different assumptions for the polarisation variables a or v. 
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Figure 1 1 : Polarisation parameters a and v measured in the helicity frame for the kinematic 
range 60 < W 1P < 240 GeV, 0.3 < z < 0.9 and P T ^ > 1 GeV, as a function of z and P T ^. 
The measurement is compared with predictions calculated in a kx factorisation ansatz [|25l and 
with calculations in CSM (collinear factorisation) at leading E51 and next-to-leading order J9]|. 
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Figure 12: Polarisation parameters a and v in the Collins-Soper frame for the kinematic range 
60 < W JP < 240 GeV, 0.3 < z < 0.9 and P T>i> > 1 GeV, as a function of z and P T ^. The 
measurement is compared with predictions calculated in a fey factorisation ansatz [25] and with 
calculations in CSM (collinear factorisation) at leading E51 and next-to-leading order [9|. 
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Inelastic J / ip Photoproduction 

[GeV 2 ] {P 2 ^) [GeV 2 ] da JP /dP^ [nb/GeV 2 ] 



1.0 - 


2.1 


1.5 


7.75 


± 


0.82 


± 


0.70 


2.1 - 


3.5 


2.7 


4.43 


± 


0.48 


± 


0.40 


3.5 - 


5.4 


4.3 


2.55 


± 


0.28 


± 


0.23 


0.4 


7 P. 
1 .0 


P. Q 
O.O 


l.UO 


i 

± 


n i q 
U.lo 


i 

± 


n m 
U. 1U 


7.6 - 


- 10.0 


8.6 


0.677 


± 


0.084 


± 


0.061 


10.0 - 


- 13.5 


11.4 


0.391 


± 


0.048 


± 


0.035 


13.5 - 


- 20.0 


15.6 


0.156 


± 


0.020 


± 


0.014 


20.0 - 


- 26.5 


22.1 


0.0509 


± 


0.0078 


± 


0.0046 


26.5 - 


- 40.0 


30.0 


0.0175 


± 


0.0029 


± 


0.0015 


40.0 - 


- 60.0 


46.0 


0.0049 


± 


0.0012 


± 


0.0004 


60.0 - 


- 100.0 


70.0 


0.00090 


± 


0.00035 


± 


0.00008 




z 


<*> 




da-fp/dz [nb] 




0.30 - 


- 0.45 


0.375 


23.4 


± 


2.6 


± 


2.1 


0.45 - 


- 0.60 


0.525 


47.6 


± 


4.7 


± 


4.3 


0.60 - 


- 0.75 


0.675 


51.3 


± 


5.0 


± 


4.6 


0.75 - 


- 0.90 


0.825 


54.2 


± 


5.6 


± 


4.9 



Table 4: Measured differential photoproduction cross sections in the kinematic range 0.3 < z < 
0.9, Px-,v> > 1 GeV and 60 < W 7P < 240 GeV as function of the squared transverse momentum 
Fyj, and the elasticity z of the J/ip meson. The bin centre values, (PtS) an d (z), are also 
given in the table. 



Inelastic J/tp Photoproduction 



W yp [GeV] 


(W 7P ) [GeV] 






cr 7P [nb] 






60 - 


- 80 


69 


0.0269 


22.9 


± 


4.1 


± 


2.1 


80 - 


- 100 


89 


0.0192 


24.1 


± 


3.3 


± 


2.2 


100 - 


- 120 


110 


0.0145 


24.0 


± 


3.0 


± 


2.2 


120 - 


- 140 


130 


0.0112 


30.3 


± 


3.6 


± 


2.7 


140 - 


- 160 


150 


0.00891 


35.7 


± 


4.3 


± 


3.2 


160 - 


- 180 


170 


0.00716 


30.4 


± 


3.9 


± 


2.7 


180 - 


- 210 


194 


0.00832 


31.7 


± 


4.2 


± 


2.9 


210 - 


- 240 


224 


0.00621 


33.8 


± 


5.6 


± 


3.0 



Table 5: Measured photoproduction cross sections in the kinematic range Pt$ > 1 GeV and 
0.3 < z < 0.9 in bins of the photon proton centre-of-mass energy W lp . The bin centre values 
(W^p) are also given in the table. $ 7 denotes the photon flux factors [52] employed in the 
photoproduction analysis using an upper Q 2 boundary of Q 2 = 2.5 GeV 2 . For the range 60 < 
W 1P < 240 GeV a photon flux factor of $ 7 = 0.1024 is calculated. 



31 



Inelastic J/ ip Photoproduction 



Pi, [GeV 2 ] (Pi,) [GeV 2 ] d<7 7P /dP^ [nb/GeV 2 ] 



0.30 < z < 0.45 



1.0 






2.0 


1.4 


1.02 


± 


0.20 


± 


0.09 


2.0 






3.0 


2.5 


0.64 


± 


0.13 


± 


0.06 


3.0 






4.5 


3.6 


0.402 


± 


0.077 


± 


0.036 


4.5 






7.0 


5.5 


0.180 


± 


0.036 


± 


0.016 


7.0 






10.0 


8.2 


0.093 


± 


0.021 


± 


0.008 


10.0 






14.0 


11.6 


0.047 


± 


0.011 


± 


0.004 


14.0 






20.0 


16.2 


0.0210 


± 


0.0052 


± 


0.0019 


20.0 






40.0 


25.0 


0.0065 


± 


0.0018 


± 


0.0006 


40.0 






100.0 


49.0 


0.00065 


± 


0.00032 


± 


0.00006 










0.45 < 


z < 0.60 










1.0 






2.0 


1.4 


2.17 


± 


0.29 


± 


0.19 


2.0 






3.0 


2.5 


1.21 


± 


0.18 


± 


0.11 


3.0 






4.5 


3.6 


0.74 


± 


0.11 


± 


0.07 


4.5 






7.0 


5.5 


0.392 


± 


0.057 


± 


0.035 


7.0 






10.0 


8.2 


0.219 


± 


0.033 


± 


0.020 


10.0 






14.0 


11.6 


0.107 


± 


0.014 


± 


0.010 


14.0 






20.0 


16.2 


0.0497 


± 


0.0084 


± 


0.0045 


20.0 






40.0 


25.0 


0.0072 


± 


0.0015 


± 


0.0007 


40.0 






100.0 


49.0 


0.00072 


± 


0.00030 


± 


0.00007 










0.60 < 


z < 0.75 










1.0 






2.0 


1.4 


2.40 


± 


0.31 


± 


0.22 


2.0 






3.0 


2.5 


1.79 


± 


0.18 


± 


0.11 


3.0 






4.5 


3.6 


1.01 


± 


0.13 


± 


0.09 


4.5 






7.0 


5.5 


0.506 


± 


0.070 


± 


0.046 


7.0 






10.0 


8.2 


0.200 


± 


0.032 


± 


0.018 


10.0 






14.0 


11.6 


0.112 


± 


0.018 


± 


0.010 


14.0 






20.0 


16.2 


0.0413 


± 


0.0076 


± 


0.0037 


20.0 






40.0 


25.0 


0.0068 


± 


0.0014 


± 


0.0006 










75 < 


z < 90 

Aj ----- \J . l/U 










1.0 






2.0 


1.4 


2.40 


± 


0.36 


± 


0.22 


2.0 






3.0 


2.5 


1.69 


± 


0.27 


± 


0.15 


3.0 






4.5 


3.6 


0.86 


± 


0.15 


± 


0.08 


4.5 






7.0 


5.5 


0.437 


± 


0.076 


± 


0.039 


7.0 






10.0 


8.2 


0.226 


± 


0.042 


± 


0.020 


10.0 






14.0 


11.6 


0.099 


± 


0.022 


± 


0.009 


14.0 






20.0 


16.2 


0.0428 


± 


0.0098 


± 


0.0039 


20.0 






40.0 


25.0 


0.0076 


± 


0.0021 


± 


0.0007 



Table 6: Measured differential photoproduction cross sections in the kinematic range 0.3 < z < 
0.9 and 60 < W JP < 240 GeV as a function of the squared transverse momentum of the J/ip 
meson in bins of the elasticity z. The bin centre values (Pt^) are also given in the table. 
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Inelastic J/ ip Photoproduction 

z (z) dcr 7P /d^ [nb] 



1.0 < P T ^ < 2.0 GeV 



0.30 - 


- 0.45 


0.375 


14.9 


± 


2.1 


± 


1.3 


0.45 - 


- 0.60 


0.525 


28.3 


± 


3.1 


± 


2.5 


0.60 - 


- 0.75 


0.675 


31.8 


± 


3.4 


± 


2.9 


0.75 - 


- 0.90 


0.825 


33.6 


± 


4.0 


± 


3.0 






2.0 < P T ,^ 


< 3.0 GeV 










0.30 - 


- 0.45 


0.375 


5.1 


± 


0.8 


± 


0.5 


0.45 - 


- 0.60 


0.525 


11.6 


± 


1.4 


± 


1.0 


0.60 - 


- 0.75 


0.675 


14.1 


± 


1.6 


± 


1.3 


0.75 - 


- 0.90 


0.825 


13.1 


± 


1.8 


± 


1.2 






3.0 < P T>lP 


< 4.5 GeV 










0.30 - 


- 0.45 


0.375 


2.60 


± 


0.42 


± 


0.23 


0.45 - 


- 0.60 


0.525 


6.05 


± 


0.73 


± 


0.54 


0.60 - 


- 0.75 


0.675 


5.71 


± 


0.71 


± 


0.51 


0.75 - 


- 0.90 


0.825 


5.32 


± 


0.80 


± 


0.48 


P T>i> > 4.5 GeV 


0.30 - 


- 0.45 


0.375 


1.10 


± 


0.20 


± 


0.1 


0.45 - 


- 0.60 


0.525 


1.30 


± 


0.20 


± 


0.1 


0.60 - 


- 0.75 


0.675 


1.11 


± 


0.17 


± 


0.1 


0.75 - 


- 0.90 


0.825 


1.30 


± 


0.24 


± 


0.1 



Table 7: Measured differential photoproduction cross sections in the kinematic range P T $ > 
1 GeV and 60 < W rp < 240 GeV as a function of the elasticity z in bins of the transverse 
momentum of the J/ip meson. 
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Inelastic J/ tp Electroproduction 



Q 2 [GeV 2 ] (Q 2 ) [GeV 2 ] da ep /dQ 2 [pb/GeV 2 ] 



3.6 - 


6.5 


4.9 


14.98 


± 


1.97 


± 


1.27 


6.5 - 


- 12.0 


8.6 


6.33 


± 


0.75 


± 


0.54 


12.0 - 


- 20.0 


15.0 


2.11 


± 


0.33 


± 


0.18 


20.0 - 


- 40.0 


26.7 


0.74 


± 


0.12 


± 


0.06 


40.0 - 


- 100.0 


53.0 


0.141 


± 


0.029 


± 


0.012 



P^ [GeV 2 ] (P^> [GeV 2 ] da cp /dP^ [pb/GeV 2 ] 



1.0 






2.2 


1.6 


15.5 


± 2.7 


± 1.3 


2.2 






3.7 


2.9 


11.0 


± 2.1 


± 0.9 


3.7 






6.4 


4.9 


8.7 


± 1.4 


± 0.7 


6.4 






9.6 


7.8 


5.90 


± 0.92 


± 0.50 


9.6 






13.5 


11.2 


3.23 


± 0.53 


± 0.27 


13.5 






20.0 


16.0 


1.69 


± 0.27 


± 0.14 


20.0 






40.0 


25.7 


0.576 


± 0.083 


± 0.049 


40.0 






100.0 


51.0 


0.055 


± 0.012 


± 0.005 






z 




<*> 




da cp /dz [pb] 


0.30 






0.45 


0.375 


150 


± 26 


± 13 


0.45 






0.60 


0.525 


158 


± 22 


± 14 


0.60 






0.75 


0.675 


280 


± 31 


± 24 


0.75 






0.90 


0.825 


239 


± 29 


± 20 


w lp 


[GeV] 


(W JP ) [GeV] 


da ep /dW lp [pb/GeV] 


60 






80 


69 


0.89 


± 0.16 


± 0.08 


80 






100 


89 


1.03 


± 0.15 


± 0.09 


100 






120 


110 


0.77 


± 0.12 


± 0.007 


120 






140 


130 


0.75 


± 0.11 


± 0.06 


140 






160 


150 


0.71 


± 0.11 


± 0.06 


160 






180 


170 


0.55 


± 0.10 


± 0.05 


180 






210 


194 


0.42 


± 0.09 


± 0.04 


210 






240 


224 


0.30 


± 0.10 


± 0.03 



Table 8: Measured differential electroproduction cross sections in the kinematic range 3.6 < 
Q 2 < 100 GeV 2 , P^ > 1 GeV and 0.3 < z < 0.9 as function of the four momentum transfer 
Q 2 , the squared transverse momentum of the J/ip meson in the photon proton rest frame Pt%-> 
the elasticity z and the photon proton centre-of-mass energy W rp . 
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Inelastic J/ ip Electroproduction 



Pfi, [GeV 2 ] (P}%) [GeV 2 ] da cp /dP^ [nb/GeV 2 ] 



0.30 < z < 0.60 



1.0 - 
4.0 - 
9.0 - 
20.0 - 


- 4.0 2.2 5.5 ± 1.1 ± 0.5 

- 9.0 5.6 3.0 ± 0.6 ± 0.3 

- 20.0 11.3 0.89 ± 0.17 ± 0.08 

- 60.0 27.0 0.11 ± 0.02 ± 0.01 


0.60 <z< 0.75 


1.0 - 
4.0 - 
9.0 - 
20.0 - 


- 4.0 2.3 3.7 ± 0.7 ± 0.3 

- 9.0 5.7 2.7 ± 0.4 ± 0.2 

- 20.0 11.3 0.92 ± 0.15 ± 0.08 

- 60.0 27.0 0.13 ± 0.03 ± 0.01 


0.75 < z < 0.90 


1.0 - 
4.0 - 
9.0 - 
20.0 - 


- 4.0 2.3 3.3 ± 0.7 ± 0.3 

- 9.0 5.7 2.6 ± 0.5 ± 0.2 

- 20.0 11.5 0.67 ± 0.13 ± 0.06 

- 60.0 27.0 0.026 ± 0.013 ± 0.002 


z (z) da cp /dz [nb] 


1.0 < P^ < 2.0 GeV 


0.30 - 
0.45 - 
0.60 - 
0.75 - 


- 0.45 0.375 60.0 ± 17.0 ± 5.1 

- 0.60 0.525 48.0 ± 11.4 ± 4.1 

- 0.75 0.675 74.6 ± 12.8 ± 6.3 

- 0.90 0.825 66.8 ± 12.9 ± 5.7 


2.0 < P^ < 3.5 GeV 


0.30 - 
0.45 - 
0.60 - 
0.75 - 


- 0.45 0.375 62.4 ± 15.1 ± 5.3 

- 0.60 0.525 67.1 ± 13.1 ± 5.7 

- 0.75 0.675 115.3 ± 16.2 ± 9.8 

- 0.90 0.825 105.0 ± 16.7 ± 8.9 


3.5 < P^ < 10. GeV 


0.30 - 
0.45 - 
0.60 - 
0.75 - 


- 0.45 0.375 28.4 ± 6.7 ± 2.4 

- 0.60 0.525 41.9 ± 7.5 ± 3.6 

- 0.75 0.675 79.6 ± 10.6 ± 6.8 

- 0.90 0.825 58.9 ± 9.5 ± 5.0 



Table 9: Measured differential electroproduction cross sections in the kinematic range 3.6 < 
Q 2 < 100 GeV 2 , P^ > 1 GeV and 60 < W rp < 240 GeV as a function of the squared 
transverse momentum in the photon proton rest frame P^L in bins of the elasticity z and the 
elasticity z in bins of the transverse momentum in the photon proton rest frame P£ 
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Inelastic J/ ip Photoproduction 
Helicity Frame 





[GeV] 


(Ptj) [GeV] 


a 


V 


1.0 
2.0 
3.0 
4.5 






2.0 
3.0 
4.5 
10.0 


1.45 
2.46 
3.65 
6.21 


+0.54 ±8:H 
-o.i5 ±8:S 
-0.18 ±S:H 

-0.28 t° H 


+Q.25 igiig 

-0.74 ±^ 

-o.o4 ig:» 

+0.59 ±°;j» 




z 




(*> 


a 




0.30 
0.45 
0.60 
0.75 






0.45 
0.60 
0.75 
0.90 


0.375 
0.525 
0.675 
0.825 


-0.65 ±8;£ 
+0.35 

-o.i8 ±8:1? 
+0.71 ±8:3 


-o.28 18$ 
+o.4o 
+0.01 ±8:1* 
-o.io ±8:S 


Collins-Soper Frame 




[GeV] 


(P T , f ) [GeV] 


Q! 


V 


1.0 
2.0 
3.0 
4.5 






2.0 
3.0 
4.5 
10.0 


1.45 
2.46 
3.65 
6.21 


+0.25 ±g:i? 

-0.26 ±°;1 5 7 

-o.o2 ±g:|g 
+o.i9^;l 


+o.4i ±g:ig 
-o.42 ±8:1? 
-o.3i 18:1? 
+Q.09 18:12 




z 




{*) 




V 


0.30 
0.45 
0.60 
0.75 






0.45 
0.60 
0.75 
0.90 


0.375 
0.525 
0.675 
0.825 


, q 47 +0.34 
-t-U.^t/ _ 28 

-o.o2 ±g : jg 
-o.oo ±8:S 
-o.o2 


-o.i8 ±8:|g 
+Q.24 ±8:1* 
-o.i6 ±8:1 
+o.5o ±8:1 



Table 10: Measured polarisation parameters in the helicity and the Collins-Soper frame as 
function of Pt,</> and z in the kinematic range Pt,v> > 1 GeV, 60 < W 1V < 240 GeV and 
0.3 < z < 0.9. 
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